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Abstract 

Quantum cluster theories are a set of approaches for the theory of correlated and disordered lattice systems, 
which treat correlations within the cluster explicitly, and correlations at longer length scales either perturbatively or 
within a mean-field approximation. These methods become exact when the cluster size diverges, and most recover 
the corresponding (dynamical) mean-field approximation when the cluster size becomes one. Here we will review 
systematic dynamical cluster simulations of the two-dimensional Hubbard model, that display phenomena remarkably 
similar to those found in the cuprates, including antiferromagnetism, superconductivity and pseudogap behavior. We 
will then discuss results for the structure of the pairing mechanism in this model, obtained from a combination of 
dynamical cluster results and diagrammatic techniques. 
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1. Introduction 

Despite two decades of intense studies, high- 
temperature superconductivity remains a mystery 
and represents one of the most important outstand- 
ing problems in condensed matter science today. 
The two-dimensional Hubbard model (Q; 0) of the 
cuprates and related models have been studied ex- 
tensively with numerical techniques, such as exact 
diagonalization and quantum Monte Carlo (QMC), 
on finite size lattices (0). These studies indicate 
that d-wave pairing correlations develop as the 
temperature is lowered (see e.g. Ref. (0)). However, 
the numerical expense or minus sign problem have 
limited finite size techniques to lattices too small 
and temperatures too high to draw conclusive evi- 
dence for the existence of d-wave superconductivity 
in these models (Jif). 

Recently, quantum cluster theories (G) have been 



used successfully to shed new light on the physics 
of the 2D Hubbard model. In these approaches a 
finite size cluster is embedded in a self-consistent 
host, designed to represent the rest of the system. 
The Cellular Dynamical Mean Field Theory |2j) uses 
a cluster in real space with open boundary condi- 
tions, while the Dynamical Cluster Approximation 
(DCA) treats a cluster in reciprocal space. Cor- 
relations within the cluster are treated explicitly, 
while longer-ranged correlations beyond the cluster 
size are described in a mean-field. In contrast to fi- 
nite size techniques, these approaches remain a non- 
trivial and physically relevant approximation to the 
thermodynamic limit even for small cluster sizes. For 
a cluster consisting of only a single site, these tech- 
niques reduce to the dynamical mean field theory 
^ which has been used successfully to investigate 
the Mott metal to insulator transition. To describe 
d-wave pairing, however, one needs to use a cluster 
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of at least four sites, which represents the mean-field 
solution for d-wave superconducting order since it 
does not allow for pairfield phase fluctuations fioi) . 
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Fig. 1. The temperature-doping phase diagram of the 2D 
Hubbard mod el obtained with a 4-site DCA calculation. 
(Jarrell et al. till) ') 

To illustrate the power of quantum cluster theo- 
ries for even small cluster sizes, we review in Fig. 1 
the temperature-doping phase diagram of a 2D Hub- 
bard model with nearest neighbor hopping t and 
Coulomb repulsion U = 8t obtained with a 4-site 
DCA calculation (Ti'). The phase diagram is remark- 
ably similar to the universal cuprate phase diagram. 
It displays an antiferromagnetic phase near half- 
filling, a d-wave superconducting phase at finite dop- 
ing and, as shown in Fig. 2, pseudogap behavior 
in the density of states and spin susceptibiliy at 
higher temperatures. Very similar results obtained 
with CDMFT (12) and cluster perturbation theory 
(a quantum cluster theor y w here the host is not 
treated self-consistently) 1)131) have confirmed this 
picture. 

In this paper, we will review DCA simulations for 
a 2D Hubbard model with nearest-neighbor hopping 
t and an on-site Coulomb interaction U. We will 
review evidence that superconductivity persists in 
larger clusters (T3) , as well as recent studies of the 
nature of the pairing interaction mediating d-wave 
superconductivity (fl^fish in this model. In the fol- 
lowing section 2 we review quantum cluster theories 
with a focus on the dynamical cluster approxima- 
tion. In Sec. 3 we explore the cluster size and shape 
dependence of the pairing correlations and the su- 
perconducting transition temperature Tc- Then, in 
Sec. 4 we investigate the nature of the pairing inter- 
action by examining the leading eigenvalue of the 
homogeneous particle-particle Bethe-Salpeter equa- 
tion and the momentum and frequency dependence 




Fig. 2. DOS for various dopings 5 = 1 — (n) in the 2D 
Hubbard model at T = 0.092t and U = 8t calculated with 
DCA/QMC for a 4-site cluster, = 4. Inset: Uniform 
spin-susceptibility as a functi on o f temperature. Energies are 
in units of 4t. (Jarrell et al. 

of the corresponding eigenfunction. Sec. 5 contains 
our conclusions. 



2. Quantum Cluster Theories: Dynamical 
Cluster Approximation 

Quantum cluster approaches share the general 
idea to approximate the effects of correlations in the 
bulk lattice with those on a finite size quantum clus- 
ter (@) . This enables a mapping of the bulk lattice 
problem to an effective cluster embedded in a self- 
consistent bath designed to represent the remaining 
degrees of freedom. In contrast to finite system sim- 
ulations which determine the exact state of a finite 
system, quantum cluster techniques give approx- 
imate results to the bulk thermodynamic result. 
In this approximation, short-ranged correlations 
within the cluster are treated explicitly, while the 
longer-ranged physics is described on the mean-field 
level. By increasing the cluster size, the quantum 
cluster theories systematically interpolate between 
the single-site dynamical mean-field result and the 
exact result, while remaining an approximation to 
the thermodynamic limit for finite cluster size. In 
the following we will focus on one of the quantum 
cluster methods, the DCA, which has been used to 
obtain the results presented in this manuscript. 

The essential assumption of the DCA is that 
short-range quantities, such as the single-particle 
self-energy S, and its functional derivatives, the 
two-particle irreducible vertex functions, are well 
represented as diagrams constructed from a coarse- 
grained propagator G. To define G, the Brillouin 
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Fig. 3. In the DCA the Brillouin zone is divided into A^c 
cells each represented by a cluster momentum K. Irre- 
ducible quantities such as the single-particle self-energy S 
and two-particle irreducible vertices are constructed from 
coarse-grained propagators G'(K) that are averaged over the 
momenta k within the cell represented by K. 

zone in two dimensions is divided into Nc = 
cells of size 2tt/L'^. As illustrated in Fig. 3, each 
cell is represented by the cluster momentum K in 
its center. The coarse-grained Green function G(K) 
is then obtained from an average over the N/Nc 
wave- vectors k within the cell surrounding K, 

G(K,c^„) = 



N ^ iujn - ej^^y. + fi- Sc(K, cj„) ' 

Here the self-energy for the bulk lattice I](K-|-k, ujn) 
has been approximated by the cluster self-energy 
Ec(K,Ci;„). Consequently, the compact Feynman 
diagrams constructed from (^(KjWn) collapse onto 
those of an effective cluster problem embedded in 
a host which accounts for the fluctuations arising 
from the hybridization between the cluster and 
the rest of the system. The non-interacting part of 
the effective cluster action is then defined by the 
cluster-excluded inverse Green's function 

g-i(K,c^„) =G-i(K,c^„) + S,(K,c^„) (2) 

which accounts for the hybridization between the 
cluster and the host. Given Q^^{K.,ujn) and the in- 
teraction on the cluster ?7 "-iT^-ji' one can then 
set up a Hirsch-Fye quantum Monte Carlo algo- 
rithm ifl^ to calculate the cluster Green's func- 
tion Gc(K,w„) and from it the cluster self-energy 
Sc(K,ct;„) which is used in Eq. (1) to re-calculate 
the coarse-grained Green function G(K,t^„) (dli). 
This process is then iterated to convergence. 

Since a determinantal Monte Carlo method is 
used to solve the effective cluster problem, there is 



also a sign problem for the doped Hubbard model. 
However, the coupling of the cluster to the self- 
consistent host significantly reduce the sign prob lem 
so that lower temperatures can be reached (|l7l ) . 

With the Hirsch-Fye QMC algortihm, one may 
also calculate the cluster two-particle Green's func- 
tion Gc2(i^4, K3;K2, Ki) with K = (K, itj„, a). Us- 
ing Gc{K) and Gc2{K4, K^; K2, Ki), one can extract 
the cluster four-point vertex F from 

Gc2iK4,K3;K2,Ki) = 

— Gc{Kl) Gc{K2) [Ski,K4Sk2,K3 — Ski,K3Sk2,K4] 

T 

+ ^ 5Ki+K2.K3+KiGc{Ki) GciKs) 

X F {Ki, K3;K2, A'l) G,(X2) Gc(ifi) . (3) 

Then, using Gc and F, one can determine the 
corresponding irreducible particle-hole and particle- 
particle vertices on the cluster, F^''' and F^p, respec- 
tively, from the corresponding Bethe-Salpeter equa- 
tions. Two-particle Green's functions for the original 
bulk lattice, G2{ki, k^, k2, ki), with k = (k, icun, c), 
are then obtained by using the cluster irreducible 
vertices, F^'' and F^p, as approximations for the lat- 
tice irreducible vertices. Lattice susceptibilities such 
as the pair-field susceptibility discussed in Sec. 3 
may then be obtained from G2(fe4, k3,k2,ki). n 

3. Superconductivity in the 2D Hubbard 
model 

With increasing cluster size, quantum cluster 
theories including the DCA progressively include 
longer-ranged fluctuations while retaining some 
mean- field character {.6). With respect to the 4- 
site phase diagram shown in Fig. 1, larger clusters 
are thus expected to systematically drive the Neel 
temperature to zero and hence recover the Mermin- 
Wagner theorem. In contrast, superconductivity 
may persist as Kosterlitz-Thouless order, even in 
the large cluster limit ifiol) . 

A measure of the strength of the d-wave pairing 
correlations is given by the d-wave pairfield suscep- 
tiblity 

P, ^ r dT{A,{r)Alm (4) 
Jo 

with AI = l/'2VNJ2i,si-'^y(^h4+si- Here, S sums 
over the four near-neighbor sites of I. In the DCA, 
Pd is calculated from the two-particle lattice Green's 
function G2(A:4, A:4, /c2, ^i) discussed in Sec. 2. The 



3 



divergence of Pd indicates a phase transition to a 
d-wave superconducting state. 
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Fig. 4. Different cluster sizes and shapes used in the DCA 
study of pairing in the 2D Hubbard model. (Maier et al. t lO)) 

To select different cluster sizes and shapes, it is 
useful to adopt the cluster selection criteria origi- 
nally introduced by Betts et al. (18) in a numerical 
study of the 2D Heisenberg model. For the Heisen- 
berg model, Betts et al. developed a grading scheme 
to determine which clusters should be used. An im- 
portant selection criterium was the completeness of 
near-neighbor shells compared to the infinite lattice. 
They found that a finite size scaling analysis was 
greatly improved when only the most perfect clus- 
ters were used. In Ref. ( 10), some of us generalized 
Belt's arguments to generate clusters appropriate 
to study superconductivity. For a d-wave order pa- 
rameter, one needs to take into account its non-local 
4-site plaquette structure. 

Fig. 4 shows the arrangement of independent d- 
wave plaquettes in the clusters. Denoting the num- 
ber of independent near- neighbor plaquettes by Zd, 
the inifinite lattice has Zd = 4. The 4-site cluster 
encloses just one d-wave plaquette, so Zd = 0. In 
this case, the effective action does not contain any 



Table 1 

Number of independent neighboring d-wave plaquettes Z^i 
and the values of T,}'" obtained from linear fits of the pair- 
field susceptibility in Fig. 5. (Maier et al. ((lOj)) 



Cluster 


Zd 




4 


(MF) 0.056 


8A 


1 


-0.006 


18A 


1 


-0.022 


12A 


2 


0.016 


16B 


2 


0.015 


16A 


3 


0.025±0.002 


20A 


4 


0.022 


24A 


4 


0.020 


26A 


4 


0.023 



pairfield fiuctuations. Thus, the 4-site cluster rep- 
resents the mean-field result for d-wave pairing and 
hence overstimates Tc- In the 8 A cluster, there is 
room for one more d-wave pair {Zd = 1). The same 
neighboring plaquette is adjacent to its partner on 
all four sides, phase fluctuations are replicated and 
hence overestimated, and thus Tc is suppressed. For 
the 16B cluster, Zd = 2, while for the oblique 16A 
cluster, Zd = Thus, one expects the pairing corre- 
lations in the 16B cluster to be suppressed relative 
to those in the 16A cluster. The number Zd of inde- 
pendent neighboring d-wave plaquettes for the clus- 
ters in Fig. 4 are hstcd in Tabic 1. The 20A, 24A and 
26A clusters all have = 4 and thus are expected 
to show the most accurate results. 



3.0 



2.5 



2.0 



9; 1.5 



1.0 



0.5 



0.0 



4A + 

8A X 
12A 
16B 
16A >^ 

20A o X 
24A • 1 

26A A X 1 ■ +■ 

x^l' 




0.02 0.04 0.06 



0.1 



0.2 0.3 
T/t 



0.5 



Fig. 5. The inverse pairfield susceptibility versus temperature 
of the 2D Hubbard model calculated w ith DCA for different 
cluster sizes and shapes. (Maier et al. llOh ) 

Fig. 5 shows results for the inverse pairfield sus- 
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ceptibility \/Pd versus temperature for U ~ At and 
(n) = 0.9. As expected the mean-field Nc = A result 
overestimates Tc, and the results for the 8 A and fSA 
clusters, both of which have Zd = 1, do not give a 
positive Tc- However, results on larger clusters with 
Zd > i fall nearly on the same curve. The results 
suggest that the 2D Hubbard model with U ~ At 
and (n) — 0.9 has a transition to a d-wave super- 
conducting phase at a finite temperature Tc- Due 
to the mean-field character of the DCA, one expect 
mean- field behavior close to Tc. The results for Tc 
of a mean-field linear fit of 1/ Pd close to the transi- 
tion are listed in Table 1. For all cluster with Zd > 

3, a transition temperature Tc w 0.023^ ± 0.002^ is 
found. 

4. The Structure of the Effective Pairing 
Interaction 

As discussed in the previous section, DCA calcu- 
lations find evidence for a finite temperature d-wave 
superconducting phase in the 2D Hubbard model. 
Here we review insight obtained by combining nu- 
merical results and diagrammatic methods to deter- 
mine the structure of the interaction mediating the 
d-wave pairing The basic idea is to focus on 

the 4-point vertex T shown in Fig. 6 calculated with a 
DCA/QMC simulation from Eq. 8. From this vertex 
and results for the single-particle Green's function, 
one can determine the irreducible particle-particle, 
r^'^', and particle-hole, F^'', vertices using the Bethe- 
Salpeter equations shown in Fig. 6a and b. Since the 
Monte Carlo results for F satisfy crossing symme- 
try, the pairing interaction F^^ may be related to the 
particle-hole channels through the equation shown 
in Fig. 6. 

In order to determine the structure of the dom- 
inant correlations in the system, the cluster ir- 
reducible vertices and the lattice single-particle 
Green's function are used to calculate the eigenval- 
ues and eigenfunctions of the homogeneous Bethe- 
Salpeter equation. For example in the particle- 
particle channel 

k' 

KMK) (5) 

with a similar equation using F'''' for the particle- 
hole channel. Here, Vpp{K\K') = rPP{K, -K- K\ -K') 
and the sum over k' denotes a sum over both mo- 
mentum k' = K' + k' and Matsubara w,,/ variables. 



r"" * r'" r 











+ 














r 


r 



Fig. 6. Bcthe-Salpeter equations for (a) the particle-particle 
and (b) the particle-hole channels showing the relation- 
ship between the full vertex V, the particle-particle irre- 
ducible vertex T^^ and the particle-hole irreducible vertex 
r'''' , respectively, (c) Decomposition of the irreducible par- 
ticle-particle vertex V^^ into a fully irreducible two-fermion 
vertex Airr plus contributions from the particle-hole chan- 
nels. All diagrams represent DCA cluster qu ant ities, includ- 
ing the Green's function legs. (Maier et al. 

By assumption, irreducible quantities like F^^' and 
(j)a do not depend on k', allowing to coarse-grain the 
Green function legs, yielding an equation that de- 
pends only on coarse-grained and cluster quantities 



K' 



withxr(^') = ^Ek'GT(K'+k',^c.„OGx(-K'- 
k', -iuJn')- 

The temperature dependence of the leading eigen- 
values in the pairing, charge density and magnetic 
channels is plotted in Fig. 7. Here, U = At and (n) = 
0.85 and the results were obtained with the quantum 
Monte Carlo dynamic cluster approximation on the 
24-site cluster shown in Fig. 4. An eigenvalue reaches 
1 when instabilities to phase transitions occur. As 
the temperature is lowered, the leading particle-hole 
eigenvalue occurs in the magnetic channel and has a 
center of mass momentum Q = (tt, tt) and frequency 
= 0. This eigenvalue saturates at low temper- 
atures. The leading particle-particle eigenvalue is a 
spin singlet, and as shown in Fig. 10, its eigenfunc- 
tion (f>(;(K, w„) has dx-2_y2 symmetry. 

The temperature dependence of the d-wave eigen- 
value Xd{T) obtained on a 4-site cluster for (n) = 0.9 
with U/t = 4, 8 and 12 is displayed in Fig. 8. The 
d-wave eigenvalue is largest when U is of order the 
bandwidth W = 8t, consistent with the notion that 
it is important to have strong antiferromagnetic cor- 
relations. Fig. 9 shows Xd{T) for various band fiU- 
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Fig. 7. Leading eigenvalues of the Bethe-Salpeter equation 
(e.g. Eq. 6) in various channels for U/t = 4 and a site occu- 
pation (n) =0.85 calculated on the 24-site cluster shown in 
the inset. The Q = (tTjTt), Wm = 0, S = 1 magnetic eigen- 
value saturates at low temperatures. The leading eigenvalue 
in the singlet Q = (0,0), cUm = particle-particle channel 
has d^2_,, 2 sy mmetry and increases toward 1 at low tem- 



peratures IKJ) . The largest charge density eigenvalue occurs 
in the Q = (0, 0), ui m = channel and saturates at a small 
value. (Maier et al. (3)) 



Fig. 8. The d^2_y2 eigenvalue \d{T) versus T/t for U = At, 
8t and 12t and (n) = 0.90. (Maier et al. IT3) ') 

ings (n) for U = 6t calculated on a 4-site cluster. 
Arf(T) increases as the system is doped towards half- 
filling. However, at half-filling the dominant eigen- 
value occurs in the particle-hole magnetic channel 
as indicated by the dashed line, and the ground- 
state has long-range antiferromagnetic order. Fur- 
thermore, for large U , a Mott gap opens leaving no 
holes to pair. 

We now turn to the momentum and frequency de- 
pendence of the eigenfunction 4>rf(K, tUn) that corre- 
sponds to the d-wave eigenvalue Xd{T). Fig. 10 illus- 
trates the momentum dependence of ^^(K, w„) for 
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Fig. 9. The d^2_y2 eigenvalue Xd{T) versus T/t for various 
band fillings (n) for U/t = 6 and Nc = 4. The dashed line 
represents the leading eigenvalue in the Q = (t t, tt ), 

S = 1 particle-hole channel at half-filling. (Maier et al. il3b ) 

ujn = ttT calculated on a 24-site cluster. Here, U = 
8t and T = 0.22i, for which \d{T) = 0.42. The val- 
ues of K lay along the dashed line shown in the inset 
of Fig. 7. Here, the d-wave structure of <I>rf(K,ct;„) 
is apparent. As shown in the inset, $d(K,Ci;„) falls 
off rapidly for momenta K moving away from the 
Fermi surface towards the zone center. 



(TT.O) 




(0,7r) 



Fig. 10. The d^2_y2 eigenvector <l?tj(K,a;n) at u)n = ttT, 
normalized to its value at K = (tt, 0), versus K for U/t = 8, 
band filling (n) = 0.9 and T/t = 0.22. In the main figure, 
the K points move along the dashed line shown in Fig. 7. 
The inset shows the b ehav ior of when K varies along the 
kx axis. (Maier et al. 0)) 

The frequency dependence of 4>d(K,w„) at the 
antinodal point K = (tt, 0) is plotted in Fig. 11 
for U/t = 4 and 8 calculated on a 4-site cluster 
at (n) = 0.9. Also shown is the spin susceptibility 
Xs(Q,a;,„) for Q = (tt, tt) normalized to coincide 
with $d(K, LUn) at w„ = ttT. One sees that the pair- 
ing interaction is retarded and characterized by the 
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same frequency scale as the antiferromagnetic spin 
susceptibility. As U increases, both $c;(K,a;„) and 
Xs(Q, (jJm) fall off more rapidly reflecting the reduc- 
tion in frequency scale set by the exchange energy. 
From these results it is clear that the dynamics of 
the pairing interaction is associated with the spin 
fluctuation spectrum. 
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Fig. 11. The Matsubara frequency dependence of 
^a(K,oj„)/<S>a(K,nT) with K = (tt, 0) for U/t = 4 and 8 
calculated for Afc=4, T/t = 0.125 and band filling (n> = 0.90. 
Also shown is the frequency dependence of the normalized 
spin susceptibility 2\s(Q,oJm) /\ys(Q, 0) +Xs(Qi27rT)] for 
Q = (vr, tt). (Maier et al. ITsI) ') 

More can be learned about the mechanism re- 
sponsible for dx2_y2 pairing in the doped Hubbard 
model, by decomposing the pairing interaction Fpp 
as shown in Fig. Ic. Here, the irreducible particle- 
particle vertex is given as a combination of a fully 
irreducible two-fermion vertex Ain- and partially re- 
ducible particle-hole exchange contributions (19). 
For the even frequency, even momentum part of the 
irreducible particle-particle vertex V^^^^{K,K') = 
1/2[TPP{K,K') + TPP{K,-K% one obtains 

+ U,iK,K') + ^^^{K,K') (7) 

with K = (K, iun)- The subscripts d and m denote 
the charge density {S = 0) and magnetic (5 ~ 1) 
particle-hole channels 

K') = 1 [Frf/„(i^ - K'; K\ ^K) 



^7/JK-K'-K',-K) 



(8) 



The center of mass and relative wave vectors and 
frequencies in these channels are labeled by the first, 
second and third arguments respectively. Using the 
Monte Carlo results for T^p, $d and $„, one can 
determine the fully irreducible vertex Ajn.. 

Fig. 12a shows the irreducible particle-particle 
vertex r^'^'(K|K') versus momentum transfer q = 
K — K'. Here, K' is set to (7r,0), K takes values 
along the dashed line in Fig. 7 and ujn — i^n' — ttT. 
When the temperature is lowered, F^'^' increases with 
momentum transfer q, as one expects for a d-wave 
pairing interaction. As can be seen from Figs. 12b-d, 
the dominant contribution to the pairing interaction 
clearly comes from the magnetic exchange chan- 
nel, <i>m, while the charge channel, ^d, and the fully 
irreducible contribution Aiir are rather independent 
of temperature and momentum transfer. 



5. Conclusions 

Quantum cluster theories treat correlations 
within a finite cluster explicitly, while approximat- 
ing those beyond the cluster at the mean-field level. 
The 4-site DCA phase diagram for the 2D Hub- 
bard model contains antiferromagnetism, d-wave 
superconductivity and pseudogap behavior in good 
agreement with experiments on the cuprates. In 
large clusters, superconductivity persists and the 
transition temperature is independent of the cluster 
size. Results for the pairing interaction F^'^'(K|K') 
show that its eigenfunction ^dCKjWn) of the lead- 
ing low temperature eigenvalue is an even fre- 
quency singlet with d^2_y2 symmetry. This implies 
that rP''(K|K') is short-ranged and increases as 
the momentum transfer q = K — K' increases. 
The frequency dependence of $d(K,Ci;„) reflects 
the frequency dependence of the spin susceptibility 
Xs(Q,a;,„) for Q — (7r,7r). The leading eigenvalue 
Xd is largest for U of order the bandwidth and in- 
creases towards half-filling, but for large U a Mott- 
Hubbard gap opens leaving no holes to pair. Finally, 
an exact decomposition of F^p shows that the dom- 
inant contribution to the pairing interaction comes 
from the 5 = 1 particle-hole channel. 
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Fig. 12. (a) The irreducible particle-particle vertex rPP(K|K') versus q = K — K' for various temperatures with u!„ = uj^i = ttT. 
Here, K = (tt, 0) and K' moves along dashed line in the 24-site cluster shown in the inset of Fig. 7. (b) The q-dependence of 
the fully irreducible two-fermion vertex Ai^r- (c) The q-dependence of the charge density (S = 0) channel for the same 

set of temperatures, (d) The q-dependence of the magnetic (S = 1) channel §3>m- (Maier et al 
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